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LONG-TERM  GOALS 

The  long-term  goal  of  this  effort  is  to  help  improve  turbulent  mixing  parameterization  in  3-D 
numerical  ocean  circulation  models  used  for  studying  the  oceans,  and  in  operational  centers,  for 
nowcasting/forecasting  the  oceanic  state. 

OBJECTIVES 

The  principal  objective  of  this  research  is  to  help  improve  second  moment  closure  (SCM)  based  ocean 
mixed  layer  (OML)  models  that  are  in  current  (and  potential  future)  use  in  Navy  community  and 
operational  ocean  circulation  models. 

APPROACH 

Extensive  research  over  the  past  three  decades  has  established  second  moment  closure  as  a  reasonable 
compromise  between  resource-intensive  techniques  such  as  large  eddy  simulations  (LES)  and  simple 
bulk  mixed  layer  models  (for  example,  Large  et  al.,  1994).  The  SMC  approach  in  its  most  practical 
form  reduces  to  a  two-equation  model  of  turbulence,  with  prognostic  equations  for  the  turbulent  kinetic 
energy  (TKE)  and  the  turbulence  length  scale  (TLS),  and  algebraic  expressions  for  the  mixing 
coefficients  {Mellor  and  Yamada,  1982;  Galperin  et  al.,  1988;  Kantha  and  Clayson  1994,  2000).  These 
so-called  algebraic  stress  closure  models  have  become  the  mainstay  of  the  US  Navy  operational  ocean 
and  atmosphere  forecast  models,  for  example  the  Shallow  Water  Analysis  and  Forecast  System 
(SWAFS)  run  routinely  at  NAVOCEANO  and  Coupled  Ocean  Atmosphere  Mesoscale  Prediction 
System  (COAMPS)  run  at  FNMOC,  as  well  as  many  civilian  operational  (NOAA  NCEP)  and  research 
(NCAR  WRF)  forecast  systems. 

However,  three  decades  of  research  and  over  a  decade  of  operational  use  have  exposed  some 
shortcomings  of  the  current  SMC-based  OML  models.  For  example,  the  popular  Mellor-Yamada  (MY) 
OML  models  in  Navy  operational  use,  have  a  tendency  to  under  predict  mixing  and  hence  overestimate 
upper  layer  currents  and  SST.  The  most  glaring  conceptual  weakness  is  the  one  related  to  the 
prescription  of  the  turbulence  length  scale.  MY  models  use  an  ad-hoc  wall  correction  to  their  TLS 
equation  (Mellor  and  Yamada  1982),  whereas  the  k-^(TKE  and  its  dissipation  rate)  model  used 
extensively  by  the  European  community  (for  example,  Rodi,  1987)  exhibits  disturbing  singular 
behavior  in  parts  of  the  parameter  space  ( Burchard  and  De leers nijder,  2001).  Another  drawback  is  the 
local  nature  of  the  closure  that  does  not  work  well  under  free  convection  conditions.  Yet  another  one  is 
ignoring  the  very  important  influence  of  surface  gravity  waves  on  mixing  in  the  upper  ocean.  None  of 
the  Navy  community  ocean  models  such  as  ROMS/TOMS,  NCOM  and  HYCOM  incorporate 
completely  surface  wave  effects;  neither  do  they  account  for  non-local  effects  under  convection. 
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Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 
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Observational  data  to  compare  with  turbulence  models  are  scarce.  Microstructure  measurements  have 
not  become  a  routine  staple  of  oceanographic  measurements  as  CTD  casts  have  been  for  decades.  This 
has  led  us  to  make  microstructure  measurements  during  NURC/NRL  2006  DART  cruises  in  the 
Adriatic  Sea,  and  LASIE  2007  (Ligurian  Air-Sea  Interaction  Experiment)  in  the  Ligurian  Sea.  We  have 
taken  part  in  the  DART  06A  and  06B  cruises  in  March  and  August  of  2006,  and  been  Scientist  in 
Charge  for  the  LASIE  2007  experiment  on  board  the  Urania  R/V,  collecting  turbulence  data  using  a 
microstructure  profder.  See  Prandke  (2005)  and  Prandke  et  al.  (2000)  for  details  of  the  microstructure 
profder  used. 

WORK  COMPLETED 

We  have  proceeded  in  the  analysis  of  acquired  microstructure  data  collected  during  the  DART  cruises 
(see  Carniel  et  al.  2006  for  details)  and  made  comparisons  with  modeled  dissipation  rate  and 
diffusivity  from  ROMS/TOMS  model.  Moreover,  new  turbulence  data  have  been  acquired  during  the 
LASIE  2007  Cruise  in  the  Ligurian  Sea,  where  the  PO  has  coordinated  a  joint  effort  together  with  two 
other  R/V  acting  as  Scientist  in  Charge. 

Lirst  comparisons  between  acquired  data  and  model  results  have  been  obtained.  Ligure  1  shows  an 
example  of  further  analysis  of  the  measurements  made  during  the  March  2006  DART06A  cruise. 
Layered  density  structure  most  likely  due  to  double-diffusive  convection  resulting  from  cold  fresh 
water  masses  over  warm  salty  ones  can  be  seen  in  the  buoyancy  frequency  and  %  profiles. 
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Figure  1.  Staircase  structures  in  the  water  column  at  Secdon-3  during  DARTO 6- A  Cruise.  The 
water  column  depth  is  27  m.  Note  the  enhanced  dissipation  rate  of  thermal  variance  at  the  interface 
between  well-mixed  layers.  Measurements  were  made  on  March  20,  2006.  Upper  panels:  (a) 
temperature  (°C)  and  salinity  (PSU,  red)  in  the  water  column,  (b)  temperature  variance  dissipation 
rate  %&  (K2s-1),  (c)  density  crQ  (Kg  m-3,  blue)  and  buoyancy  frequency  (s-1).  Bottom  panels:  (d)  heat 
diffusivity  K9  (m2s-l),  (e)  TKE  dissipation  rate  (W/Kg),  (f)  eddy  diffusivity  KM  (m2s-l).  Color  lines 
are  the  average,  the  two  thin  black  lines  represent  the  bootstrap  confidence  limits. 


Figure  2  shows  the  corresponding  water  masses  depicted  at  different  depths. 
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Section  3 


Figure  2.  T/S  diagrams  for  Section-3,  with  water  masses  depicted  at  different  depths . 


Water  masses  at  Section-3  are  characterized  mainly  by  a  slightly  saltier  (38.52-38.7  PSU)  and  much 
warmer  (13.45-13.85  °C)  water  mass  w.r.t.  most  of  the  water  of  the  other  measurements,  thus 
originating  a  lighter  (29.06-29. 18  kg/m3),  representing  an  off-shore,  open  ocean  water. 

For  double  diffusive  convection  with  a  temperature  and  salinity  changes  of  AT  and  AS  across  the 
interface,  and  hence  a  density  ratio  Rp  =  (a=0.216-10'3  °C'1  and  J3=0. 77-10°  PSU'1),  the  values 

of  momentum,  salt  and  thennal  diffusivities  are  (Marmorino  and  Caldwell,  1976): 
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where  v=l  .5-10'  ms'  is  the  molecular  kinematic  viscosity  of  water. 
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For  the  two  interfaces  at  station  GS2,  AT=  0.67°C  and  AS=  0.33  PSU  so  that  Rp  =  0.56  ;  according  to 

Eq.  (1),  we  have  =7.  IT  O'4,  Ks  =  KM  =  1.4-1  O'4  mV1.  These  values  are  consistent  with  the 
diffusivities  presented  in  Kelly  (1984)  and  while  the  peak  values  in  measured  Ke  are  much  higher,  the 
average  value  for  the  water  column  is  about  2- 1 0‘5  in  rough  agreement  with  the  value  cited  above. 

For  the  bottom  salt  finger  interface  of  Section-3,  AT=  0.12°C  and  AS=  0.03  PSU  so  that  Rp  =  1.12  ; 

therefore,  according  to  Eq.  2),  we  have  Ks=  KM=  9.4- 10'4,  Ke=  6.6-1 0"4  mV1.  The  values  for  Ke  and 
Km  in  Fig.  4  are  computed  assuming  conventional  turbulence  and  hence: 
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where  Rf  is  the  flux  Richardson  number  (i.e.  the  efficiency  of  conversion  of  TKE  into  potential  energy, 
~0.2)  and  N  is  the  buoyancy  frequency.  As  can  be  seen  in  Figure  1,  peak  values  are  both  roughly 
around  10~4  mV1. 


Instead,  if  we  use  the  relationship  proposed  by  McDougall  (1988)  to  compute  the  effective  vertical 
eddy  diffusivity  of  salt  in  salt  fingers: 
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where  fR  is  the  flux  ratio.  In  this  case,  the  values  for  KM  presented  in  Figure  1  would  be  considerably 

R-  i 

higher,  since  the  coefficient  of  — — —  is  around  1,  if  we  take  a  value  of  0.9  for  the  flux  ratio  at 

1 —  f R 

Rp  =1.12  (however,  note  that  there  exist  considerable  uncertainty  in  extrapolating  the  flux  ratio 
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measurements  of  Takao  and  Narusawa  (1980)  to  values  of  R  slightly  higher  than  1).  If  we  use  the 

second  relationship  in  Eq.  (4)  the  values  would  be  similar  to  those  of  Kg  but  25%  higher.  Either  way, 
the  results  are  well  within  the  confidence  limits  of  the  bootstrap  analysis  of  the  measured  values. 

RESULTS 

When  the  turbulence  is  generated  by  both  the  momentum  flux  and  a  destabilizing  buoyancy  flux,  the 
TKE  dissipation  rate  £  in  the  mixed  layer  can  be  taken  to  be  the  sum  of  the  rates  due  to  shear-driven 
and  buoyancy-driven  turbulence. 

The  following  expression  for  scaling  can  be  used: 

£  =  £c+£s  z<D 

=  £,  z>D 

(5) 


where  in  the  convective  mixed  layer  we  have 

£c  =  J b0  z  =  0 

=  039 Jh0  0.  ID  <  z  <  0.9D 

=  0  z>D  (6) 

and  in  the  wind-stress  driven  region 

£s  =  ul  /  (a z)  0  <  z  <  0.3D 

=  3.33m* /(aD)  0.3D<z<D 

=  0  z>D  (7) 

Below  the  mixed  layer  depth  D,  the  expression  £t  =  CKL2TN 3  (9),  resulting  from  assuming  that  the 
Ozmidov  scale  Lo  is  proportional  to  this  overturn  scale  LT  ,  i.e.  L0  =  c0LT,  can  be  used. 

Considerable  effort  has  been  expended  in  determining  the  value  of  Co  (See  Thorpe  2005  for  a  summary 
and  Dillon  1982  for  a  lucid  discussion  of  the  Thorpe  scale  and  related  issues).  Plots  of  the  observed  LT 
versus  L0  in  the  ocean  show  considerable  scatter  ( Osborne  1980,  Dillon  1982,  Wesson  and  Gregg 
1994,  Mourn  1996,  Peters  et  al.  1988,  Ferron  et  al.  1998,  Gargett  1999,  Caldwell  1983,  Galbraith  and 
Kelly  1996,  Fer  et  al.  2004)  typical  of  ah  microstructure  measurements,  and  so  the  value  of  cq  does 
involve  some  uncertainty.  The  scatter  may  very  well  be  due  to  the  noise  and  the  resolution  of  the 
sensors,  but  the  bin  size  used  in  determining  L  j  may  also  be  an  important  factor.  In  Wesson  and  Gregg 
(1994)  measurements,  Co  ranges  between  0.25-M.0,  yielding  a  value  of  Ck  between  0.0625^-16. 

In  order  to  further  investigate  the  validity  of  the  relationship  expressed  by  £l  =  CKL2TN 3 ,  the  values  of 

the  Thorpe  scale  (i.e.  the  RMS  of  the  Thorpe  displacements)  at  different  bins  resolution  (0.5  m,  1  m,  2 
m  and  5  m)  were  computed,  using  the  full  dataset  available  at  B90  station,  collected  during  the  DART- 
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06B  cruise,  August  2006.  Corresponding  values  for  the  proportionality  constant  cq  were  around  1., 
close  to  the  lower  range  proposed  by  Wasson  and  Gregg  (1994).  The  value  of  the  proportionality 
factor  appears  to  be  dependent  on  the  length  of  the  bin  considered  and  varying  if  computed  at  different 
depths,  representative  of  different  processes. 

Figure  3  presents  a  log-log  diagram  of  the  values  of  the  LT  vs  Lo,  computed  following  Eq.  (9),  for  the 
5  m  bin  example  (average  co=1.15). 


Figure  3.  LT  and  L0  values  for  B90  values,  5  m  bin  average. 

The  two  black  lines  represent  the  LT=4L0  (upper  one)  and  LT=1/4L0 


IMPACT/APPLICATIONS 

Accurate  depiction  of  many  quantities  of  interest  to  worldwide  naval  operations,  such  as  the  upper 
layer  temperature  and  currents,  requires  accurate  simulation  of  turbulent  mixing  in  the  water  column 
and  accurate  tidal  forcing.  Operationally,  this  contributes  to  better  counter  mine  warfare  capabilities 
through  better  and  more  accurate  tracking  of  drifting  objects  such  as  floating  mines.  Other  drifting 
materials  such  as  spilled  oil  are  also  better  tracked  and  counter  measures  made  more  effective.  Other 
applications  include  search  and  rescue.  Turbulence  data  collected  under  this  project  can  help  assess 
turbulence  parameterization  in  OML  models. 
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